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Abstract 


This paper concentrates on the application of optimal control theory to 
highlight some aspects of Turkish economy. First the setup is given for 
Turkey to grow over the balanced path. Then the optimal control 
problem is identified. The control and state variables are mentioned. The 
objective is the maximization of life-time discounted utility of the society 
through optimal choice of consumption which automatically determines 
investment. We make use of Bellman’s principle to guarantee optimality. 
We make necessary assumptions (technical assumptions) to make use of 
calculus techniques for a solution. Some functions to represent utility and 
production are specified. I used the econometric techniques to estimate 
some parameters of the functions to decide upon the optimal level of 
investment for steady-state in Turkey over the period including 2001 
crisis. The corresponding differential equations are obtained as a result of 
the Hamiltonian defined. The phase diagram is prepared to analyse 
different trajectories. 
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Introduction 


Literature on growth holds a vast volume in macroeconomic theory. This subject is at the center 
of many macroeconomic research discussions. Contemporary theory on growth can be traced 
back to Solow (1956) for his contributions to the so-called exogeneous growth theory. He has 
made use of the basics developed by even Adam Smith. Later contributions to theory formalized 
by Solow were towards repairing some shortcomings of this theory on several fronts and these 
contributions have led to the endogeneous growth where the saving rates of the consumers can 
somehow be internalized within the economy. This innovation in theory made the framework 
much closer to real consumption behaviors of individuals and households. 


Much earlier studies on growth belong to the classical school of thought. Adam Smith (1904) 
proposed that growth is supply-sided and total output is a function of labor, capital and land. 
Therefore, growth was stated to be driven by population, investment and land progress as well as 
technological advances. Smith thought that all these factors were endogenous and he had given 
emphasis to specilaization as a particular engine to fuel growth. He suggested that there is a limit 
to growth as a result of an ultimate halt of capital accumulation and population increase. 


The other major profile in the classical school of thought, David Ricardo, made a remarkable 
contribution to the discourse of Smith by saying that land is limited and cannot be cultivated as 
Smith suggested. And Malthus had stated that population growth depends on the food supply 
and would ultimately limit economic growth thereby leading to the poorness of individuals. 


The neoclassical school of thought could not add much to the pile of the classical school. It was 
early 1900s at which Schumpeter completed his Theory of Economic Development and had 
mentioned the steady state as his reference point. Schumpeter assumed population growth was 
exogenous as opposed to Smith and has taken the enterpreneurs as the main driving force of 
development (Schumpeter’s arguments were around development instead of growth). 


It was the middle of the 20th century when Robert M. Solow had presented his growth model 
which became famous by his name, Solow or neoclassical growth model. Solow appeared as the 
most prominent economist in the growth literature. He proposed a model that handled the capital 
output ratio fixed and this ratio guranteed the system to move back to steady state. He assumed 
that there are CRS, perfect competition, perfect information flow and no externalities. These 
were technical assumptions instead of being realistic. He had to make these assumptions to 
guarantee the steady state mathematically. And this is the reason of why Solow’s Model lacked 
empirical evidence in many studies. Furthermore, Solow had assumed that the individual savings 
and technology were exogeneous. In the same year Swan had made a similar study so that the 
model is sometimes called the Solow-Swan Model. Solow had introduced the concept of growth 
accounting through which growth is decomposed into its capital accumulation, labor force 
increase and the residual that handles all others including technology. 


There are two main shortcomings of the neoclassical growth models. First of all, it lacks 
empirical evidence. Many of the empirical studies have disproved the model. Secondly, 
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technology, which is the key determinant of growth was taken exogeneous. The model could no 
longer be used. Therefore, the model was modified by two economists, Kass (1965) and 
Koopmans (1967). They had made use of the earlier studies by Ramsey (1928) who had taken 
the saving rate endogenous. That is why the model is known as the Ramsey-Kass-Koopmans 
Model. 


The progress of the growth theory in last decades are mainly due to the studies of Romer (1986) 
and the Nobel Laurate Robert E. Lucas (1988) that popularized the earlier study of Romer. This 
streamline of growth theory is called “New Growth Theory” and questions the results and 
assumptions of the previous one. Basically, the returns to scale are taken increasing and the 
competition is assumed imperfect which reflect the facts of real economies in a better way. 


This paper superficially summarizes the very well-known literature of the growth theory in the 
first section, Introduction. The second section introduces the model and the basics of the optimal 
control theory to resolve this model. Third section includes the phase diagram based on the 
optimal control analysis. Section 4 estimates the parameters of the econometric model which 
handles the macroeconomy of Turkey. Section 5 lists down the concluding remarks. 


Optimal Control Problem of Growth 


We assume a production function with the arguments of capital (K), labor (L), and others, (A). 
YO) = FKA, LOAM) [1] 


We have the population growth rate of n and technology growth rate of g. Both growths are 


continuous: L(t) =e" and A(t) =e". 
What is not depreciated and not consumed out of production is the increment of the capital stock: 


K(t) = F(K(@), LAW) — 8K) — CA) [2] 


Letting c(t) = co Z(t) = ALi 


_ Kt) 
L(t) AY 


= we can write the equation in its 
A(t) L(t) 


intensive form: 


K(t) = f(K@)-Z)—-(n+6 + g)«) [3] 
The society wants to maximize the lifetime utility function where the life of the society is known 


to be infinity: 
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co 


max v) fe "U(c@)at [4] 


0 


where / is the discount factor. 


We assume that the utility function is in the form of Constant Relative Risk Aversion (CRRA): 


1-0 
Uc) =~ [5] 


when 0+1 and when @=1 we assume that U(c(t)) = Inc(t). There is no empirical research to 


figure out the numerical value of @ to the best of my knowledge. The numerical value is 
assumed to be 0.355. This number follows from the empirical research of Holt and Laury (2002). 


Substituting this utility function into the maximization problem leads to: 


max p |e” aly dt 
ae 
Substituting Z(t)e*' for c(t) and rearranging give: 
j VAO 


max [e EnD Z dt [6] 
J E 
Defining the utility function U in terms of Z , and using a to handle the power of e the problem 
turns out to be: 
max Í e“U(Z(t))dt [7] 
0 

Both output and consumption levels can be figured out when « and Z are worked and figured 
out. We have the typical social planner in charge one more time. The social planner is directly 
responsible for the life-time discounted utility maximization of the whole society. 


So that the social planner has the problem of: 
max fe “U(Z())dt 


0 
st 
Kt) = fE- ZO- n+ 6+ gE) 
and 
K(0)= Ky 
This problem is much different than an ordinary maximization problem because the function to 
be maximized is the integral and the constraint is a differential equation. This is a typical optimal 


control problem. Pontryagin is the prominent Russian mathematician to study these kinds of 
problems. Dorfman (1969) applied the theory to economic problems. One can read Kamien and 
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Schwarts (1981), or Stokey et. al. (1989) for very complete exposure of the material on optimal 
control theory. 


This problem can be solved with the help of Pontryagin’s Maximum Principle. The allocation 
that is found as a result of this problem is at the same time Pareto Optimum. The state variable is 
x(t) and the control variable is Z(t). We define the co-state variable 1 to form the present value 


Hamiltonian, H . 


H(k,Z,A,t) =e “U(Z(t) + Af (KAO -ZA — (n+ 6 + g)K(t) [8] 


The sufficient conditions for an optimal solution are: 
A(t) =e “U'Z(t) [9a] 


AQ) =—f'(KM)—-(+ 6+ ga) [9b] 


lim A(t)«(t) =0 [9c] 


In order to solve this system of differential equations one would rather get rid of the co-state 
variable, /, first. Differentiating Eq. (9a) with respect to time gives: 
A(t) =U'(Z(T)) Z(t) —ae“"'U' Z(t) [10] 
Dividing both sides of this equation by both sides of Eq. (9a): 
At) ZOU"ZO 
A(t) U'Z(t) 
Incorporating Eq. (9b), and multiplying both sides of the equation by Z(t): 


Z(t) = pif (sto) —(n+6+g+a))Z(t) [11] 


ZOU" ZO) 
UZO) 


Here @ is substituted for 


Now the two differential equations are obtained, Eq. (11) above, and Eq. (12) below. 


K(t) = f(K)-Z—(n+ 6+ g)«(t) [12] 


Phase Diagram 


Phase diagrams are very beneficial tools that are employed in analyzing the behaviors of 
functions that appear in systems of differential equations especailly in cases where the set of 
differential equations cannot be solved directly. In our situation the set of equations include some 
functions that are not identified, so we have no chance to solve these differential equations. But 
we know some characteristics of them to draw some conclusions about them over the phase 
diagrams. 
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Let’s start by writing down the two differential equations in hand one more time: 


AOE S'E- ++ +a))Z(t) [11 


K(t) = fK) -ZA -(n+8 + g)K(t) [12] 


We assign Z to vertical axis and «x to horizontal axis. Let’s indicate the steady-state values of Z 


and x by Z* and x , respectively. 


Let’s go through Eq. (11) first. Z =0 when the right-hand-side (RHS) of this equation is equal 
to 0. This is only when f'(x(t))—(n+64+g¢+a)=0 or f'(K()) =(n+6+2+a)Z(t). This is 


when x(t) is at its steady-state condition, x(t) =x". We draw a strict line at «K(t)=K to 


indicate the points at which Z (t)=0. 


What happens to Z(t) off this line? Now let’s look at Eq. (11) one more time. Let’s pick a point 


«(t)>« over the phase-diagram plane. For this point f'(x(t)) is smaller since f" is a 


decreasing function of x(t). That is why f'(«(t))—(n+6+ g +a) is negative and Z(t) (which is 
equal to f'(x(t))—(n+6+ g +a)) is also negative since we are in the Ist quadrant and Z(t) is 
positive for sure. This means Z(t) is declining (since the derivative of Z(t) is negative). That is 


why we put downward directed arrows to indicate this decline in Z(t). Similarly when 


K(t)< K”, Z(t) >0. And this time the directions of the arrows are reversed, ie. they are directed 


upward. 


Coming to Eg. (12), x®©=0 when  f(K(t)—-Z()-(n+d+g)K()=0, or 
Z(t) = f(K(t)) -—(n+ 0+ g)k(t). Z(t) =O when x«(t)=0. So that the shape of the curve for 


which «(t) =0, passes through the origin. 
Let’s think about the derivative of Z(t) = f(xk(t))—(n+6+g)x(t) to see how the function 
behaves. Since f is concave in x(t) the slope will be positive. But the slope will decline for higher 


values of x(t). So that the shape of the curve is like an unusual bell. The maximum arises at «* 


which is larger than the steady state value of x(t), « . This gives us the shape of Figure 1. 


Coming to the dynamics on some sides of this bell-shaped curve, let’s consider a point which is 
above the borderline, for this point Z(t) > f(k(t))-—(n+ 6+ g)K(t)so that 


K(t) = f (k(t) -—Z(t)-—(n+ 6+ g)x(t) <0, and x(t) declines. We have the arrows directed to 
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left. For the points surrounded by the bell-shape we have Z(t) < f(«K(t))—(n+6+ g)«(t) and 


K(t) > 0. The directions of the arrows at this region will be directed to right. 


These arrows indicate where the economy will let consumption and capital stock to move. As the 
figure indicates there is just one saddle-path that leads the economy to the steady-state. This 
saddle-path is shown by the trajectory of the additional arrows, ie from top-right to bottom-left. 
If the economy starts at any other location it will go astray. 


Figure 1. Phase diagram for Z(t) and x(t). 


If Turkey is positioned over the trajectory then she will be enjoying the travel over the saddle- 


path towards the golden route of growth. The main characteristic of the intersection of Z(t) =0 


vertical line and K(t) = Qbell-shape is that at this intersection both Z(t) and «(t) will have no 


tendency to change over time. 


In this last section of the paper we estimate parameters to intend to apply the theory to Turkish 
economy around her 2001 crisis. We start by stating the additional condition of the steady state: 


Z(t) =O and x(t) =0. 


Now let’s rewrite the essential equations under this condition: 
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Z(t) = URGO) _(n+d+g+a)Z) 
SO; 
O= f(x) -+6+¢+a) 
and; 


f' (k(t) =n+6+g+<a [lla] 
If one assumes the Cobb-Douglas type production function in its intensive from then: 


f(«K)=K*% and f'(k) =a". 
Substituting this in Eq. (11a): 


a- 


akx*'=n+d+g+a 


So 


This is the steady-state level of capital stock. 
Indeed, the consumption maximizing level of capital stock and the steady-state capital stock are 
different than each other. The consumption maximizing capital stock is the projection of the peak 


of the bell shape which is higher than the level of capital stock at the intersection of Z(t) =0 


and x(t) =0 trajectories apparently. 


Estimation of Parameters 


In this part of the study, I implement the finding to the Turkish economy over one of the severest 
crisis of its last century, namely the 2001 crisis. One has to numerically evaluate all parameters 
appearing on the RHS of Eq. (13) in order to assess the balanced growth path of Turkey over the 
mentioned crisis. The purpose of this last section is this assessment. 


Since L(t) =e”, taking logarithms of both sides: In L(t) = nt . This kind of a specification of the 
function automatically normalizes the population at time 0 to 1, ie. L(0) =e” =1. This equation 
resembles the simple regression function of type: Y = BX +e where Y corresponds to In L(t), 
Ø corresponds to n and the error term is assumed to satisfy the classical assumptions of 
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ordinary regression. But since what we are handling is time-series data we have to take care of 
stationarity. 


We gather annual data from the official website of the Central Bank, www.tcmb.gov.tr, from 
1968 to 2004 for this estimation. We form a new series by setting 1968 to 0 and that year’s 
population to 1. The series proves to be stationary as consequence of the Augmented Dickey- 
Fuller Unit Root Test executed by EViews. The estimated value of n appears to be 0.02222. The 
Eviews outputs are postponed to Appendix 1. 

The depreciation coefficient is assumed to be 0.05 since we assume that the capital used in 
manufacturing depreciates completely in 20 year’s time. 


We estimate g just like how we estimate n. Recall, A(t) =e“. After some similar arrangements 
and manipulations, we estimate Y = X +€, one more time for this case. We use the National 
Technology Index, as a proxy for A(t). The original data set is obtained from the same web site 


and raw data set is offered on a daily basis over 12.11.2001-13.9.2005. We make the necessary 
arrangements to have the consistency of data sets about their frequency. Again the series proves 
to be stationary as a result of the Augmented Dickey-Fuller Test. The regression output reveals 
that gis estimated to be 0.000232 and is highly significant, ie. p-value is 0.000. The EViews 
outputs are in Appendix 2. 


Coming to the estimation of a that stands for p — g(1—@). Since we assumed that p = 0.90 and 
@=0.355 and estimated 2 = 0.000232 , the estimate of a turns out to be 0.90 due to the very 
small numerical value of p : 


The last parameter to be estimated is @ of the Cobb-Douglas Production Function in its intensive 
form. 


Remembering: 
Y = F(K,L) = K" (AL)™ is written in its intensive form as: y = «* where K = ae The rest is 


arranging data and running a similar regression of In yover In«. The coefficient found is the 


estimate of a . Both series proved to be stationary at their levels. OLS estimate of œ turned out 
to be approximately 0.98 where t-statistic showed that the coefficient is highly significant. The 
regression results are given in Appendix 3. Table 1 below lists down numerical values of 
parameters to estimate « in Eq. (13). 


Table 1. Numerical values of parameters, by assumption or estimation. 


n ô g a a 


0.0222 | 0.0500 | 0.000232 | 0.90 | 0.98 


Substituting all these numerical values give the optimal level of capital per ideal labor force, x”. 


18 


Orhan, M. / Journal of Management, Economics, and Industrial Organization, Vol.1 No.1, 2017, pp.10-23. 


Concluding Remarks 


e This paper attempted to apply the optimal control theory to figure out the level of 
consumption, and thereby investment, to the economy of Turkey around her 2001 crisis. 
Real data sets of Turkey are used to find the optimizing policies. 

e Unfortunately, there are not many theoretical applications of macroeconomics to Turkish 
economy. One has to recall significant contributions of Prof. Merih Celasun (2002) as 
well as Metin-Ozcan, et. al (2001). 

e One can make more explicit analysis after completely specifying the utility function, 
production function and all parameters belonging to both of these functions. The growth 
accounting can then be handled to observe the contributions of factors of production to 
growth individually. This issue stays as an open area for further research. 

e In the last section of the paper some parameters are estimated and some are assumed to 
have certain numeric values in order to conclude « which leads to the conclusion of Z*, 
consumption. These estimates can be used by policy-makers in various developments of 


policies. 
e Last but not the least, theoretical findings of the paper can be applied to all developing 
countries. 
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Appendix 1. Unit Root Test results and the regression output: 


ADF Test Statistic -2.146026 1% Critical Value* -3.6289 
5% Critical Value -2.9472 
10% Critical Value -2.6118 


*MacKinnon critical values for rejection of hypothesis of a unit root. 


Augmented Dickey-Fuller Test Equation 
Dependent Variable: DILABOR) 
Method: Least Squares 
Sample(adjusted): 3 37 


Included observations: 35 after adjusting endpoints 


Variable Coefficient Std. Error t-Statistic Prob. 
LABOR(-1) -0.012265 0.005715 -2.146026 0.0396 
D(LABOR(-1)) 0.032330 0.174769 0.184985 0.8544 

C 0.025305 0.005154 4.909369 0.0000 

R-squared 0.148377 Mean dependent var 0.020984 
Adjusted R-squared 0.095151 S.D. dependent var 0.007146 
S.E. of regression 0.006797 Akaike info criterion -7.062774 
Sum squared resid 0.001478 Schwarz criterion -6.929459 
Log likelihood 126.5985 F-statistic 2.787656 


Durbin-Watson stat 1.999600 Prob(F-statistic) 0.076552 


Dependent Variable: LABOR 
Method: Least Squares 
Sample: 1 37 

Included observations: 37 


LABOR=C(2)*TIME 


Coefficient Std. Error t-Statistic Prob. 

C(2) 0.022221 0.000168 132.2684 0.0000 
R-squared 0.991312 Mean dependent var 0.406809 
Adjusted R-squared 0.991312 S.D. dependent var 0.229443 
S.E. of regression 0.021387 Akaike info criterion -4.825444 
Sum squared resid 0.016466 Schwarz criterion -4.781906 


Log likelihood 90.27071 Durbin-Watson stat 0.109222 
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Appendix 2. Unit Root Test results and the regression output: 


ADF Test Statistic -1.845245 


1% Critical Value* -3.4399 
5% Critical Value -2.8650 
10% Critical Value -2.5686 


*MacKinnon critical values for rejection of hypothesis of a unit root. 


Augmented Dickey-Fuller Test Equation 


Dependent Variable: D(TEK) 
Method: Least Squares 
Sample(adjusted): 6 999 


Included observations: 963 


Excluded observations: 31 after adjusting endpoints 


Variable Coefficient 
TEK(-1) -0.009073 
D(TEK(-1)) -0.037789 
D(TEK(-2)) 0.022148 
D(TEK(-3)) 0.045252 
D(TEK(-4)) -0.026804 

C -0.005615 

R-squared 0.008484 
Adjusted R-squared 0.003304 
S.E. of regression 0.024455 
Sum squared resid 0.572322 
Log likelihood 2210.196 
Durbin-Watson stat 1.997377 


Dependent Variable: TEK 
Method: Least Squares 
Sample: 1 999 

Included observations: 973 
Excluded observations: 26 


TEK=C(1)+C(2)*TIME 


Coefficient 

C1) -0.757620 

C(2) 0.000232 
R-squared 0.164580 
Adjusted R-squared 0.163719 
S.E. of regression 0.149744 
Sum squared resid 21.77310 
Log likelihood 467.9312 


Std. Error t-Statistic Prob. 


0.004917 -1.845245 0.0653 
0.032409 -1.165991 0.2439 
0.032326 0.685159 0.4934 
0.032131 1.408356 0.1594 
0.032124 -0.834412 0.4043 
0.003241 -1.732292 0.0835 


Mean dependent var 0.000183 
S.D. dependent var 0.024495 
Akaike info criterion -4.577769 
Schwarz criterion -4.547426 
F-statistic 1.637783 
Prob(F-statistic) 0.147277 
Std. Error t-Statistic Prob. 
0.009768 -77.55940 0.0000 
1.67E-05 13.83073 0.0000 
Mean dependent var -0.639958 
S.D. dependent var 0.163747 
Akaike info criterion -0.957721 
Schwarz criterion -0.947689 
Durbin-Watson stat 0.027124 
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Appendix 3. Unit Root Test results and the regression output. 


ADF Test Statistic -2.922172 1% Critical Value* -4.0113 
5% Critical Value -3.1003 
10% Critical Value -2.6927 


*MacKinnon critical values for rejection of hypothesis of a unit root. 
Augmented Dickey-Fuller Test Equation 

Dependent Variable: D(LNY) 

Method: Least Squares 

Sample(adjusted): 3 16 


Included observations: 14 after adjusting endpoints 


Variable Coefficient Std. Error t-Statistic Prob. 
LNY(-1) - 1.340020 0.458570 -2.922172 0.0139 
D(LNY(-1)) 0.174648 0.302724 0.576922 0.5756 

C 28.01735 9.582192 2.923898 0.0138 

R-squared 0.577508 Mean dependent var 0.017310 
Adjusted R-squared 0.500691 S.D. dependent var 1.283887 
S.E. of regression 0.907218 Akaike info criterion 2.830542 
Sum squared resid 9.053495 Schwarz criterion 2.967483 
Log likelihood -16.81379 F-statistic 7.517984 
Durbin-Watson stat 2.065442 Prob(F-statistic) 0.008750 
ADF Test Statistic -2.579003 1% Critical Value* -4.0113 
5% Critical Value -3.1003 

10% Critical Value -2.6927 


*MacKinnon critical values for rejection of hypothesis of a unit root. 
Augmented Dickey-Fuller Test Equation 

Dependent Variable: D(LNK) 

Method: Least Squares 

Sample(adjusted): 3 16 


Included observations: 14 after adjusting endpoints 


Variable Coefficient Std. Error t-Statistic Prob. 
LNK(-1) -0.997358 0.386722 -2.579003 0.0256 
D(LNK(-1)) 0.163368 0.289323 0.564658 0.5836 

C 21.27366 8.252885 2.577724 0.0257 

R-squared 0.445479 Mean dependent var 0.033733 
Adjusted R-squared 0.344657 S.D. dependent var 1.671197 


S.E. of regression 1.352888 Akaike info criterion 3.629770 
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Sum squared resid 20.13337 
Log likelihood -22.40839 
Durbin-Watson stat 1.943716 


Dependent Variable: LNY 
Method: Least Squares 
Sample: 1 16 

Included observations: 16 


LNY=C(2)*LNK 


Coefficient 

C(2) 0.977570 
R-squared -2.477597 
Adjusted R-squared -2.477597 
S.E. of regression 1.489377 
Sum squared resid 33.27365 
Log likelihood -28.56043 


Schwarz criterion 
F-statistic 


Prob(F-statistic) 


Std. Error t-Statistic 
0.017414 56.13713 


Mean dependent var 
S.D. dependent var 
Akaike info criterion 
Schwarz criterion 


Durbin-Watson stat 


3.766711 
4.418471 
0.039043 


Prob. 


0.0000 


20.93775 
0.798666 
3.695054 
3.743341 
1.774899 
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